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Abstract 

The MSSM predictions at the one-loop level for the weak-magnetic dipole mo- 
ments of the t lepton and the b quark are analysed. The entire supersymmetric 
parameter space is scanned with the usual GUT constraint and common squark 
and slepton mass parameters. The real part of a^f is dominated by the chargino 
contribution, being the same order as the SM one or even larger in the high tan (5 
region whereas the imaginary part, due to Higgs boson diagrams, is negligible com- 
pared to the SM value. The real part of a^ 7 is controlled mainly by charginos and 
also by gluinos, when the mixing in the bottom squark sector is large, to yield, 
for high tan/3, a contribution one order of magnitude larger than the pure elec- 
troweak SM value but a factor five smaller than the standard QCD contribution. 
The imaginary part of is the same order as in the SM. 
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The investigation of the electric and magnetic dipole moments provides very accurate 
tests of the quantum structure of the Standard Model (SM) (see |[] and references therein) 
and of possible extensions (see |3[ and references therein). The fermion- Z boson vertex, 
in higher order, involves coupling terms which in their Lorentz structure are analogous 
to the magnetic and electric dipole terms of the fermion-photon vertex with valuable 
informations on the CP-conserving and CP-violating sectors of the SM. In addition, the 
magnetic weak dipole terms exhibit non-decoupling properties 0] when the top quark is 
involved in the quantum contributions. Together with the chirality-flipping character of 
the dipole form factors one thus might expect also some insight into the mechanism of mass 
generation. Extensions of the minimal model by new renormalizable interactions influence 
the dipole form factors by non-standard loop contributions. In this article we consider the 
Minimal Supersymmetric Standard Model (MSSM) as a particularly interesting extension 
of the SM and investigate the impact of the MSSM one-loop contributions to the weak- 
magnetic dipole moments of heavy fermions. 

The most general Lorentz structure of the vertex function that couples a Z boson and 
two on-shell fermions (with outgoing momenta q and q) can be written as 



If// = ie 7jU 



7 y _ V f )-( Fa - ^ W 
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+ (? - Q)»[Fm + F El5 ] -(q + g) M [F s + F p75 ]| , (1) 

where Vf = (l( — 2s 2 v Qf), a f = H (the weak isospin of the fermion) and sfy = 1 — dyy = 
sin 2 6w. The form factors Fi(s) are defined splitting off the tree level SM terms, and 
depend only on s = p 2 = (q + q) 2 . Conservation of the vector current (U(l) gauge 
invariance) constrains the scalar form factor F$ to be zero at s = M§. The axial form 
factor Fp(M§) must also vanish for massless fermions. The vector and axial- vector form 
factor, Fy and F&, are the only ones that have to be renormalized while the others are 
finite. The form factors Fu and Fe are related to the weak dipole moments of the fermion 
/ with mass m/ as follows ||: 

AWMDM = aj = -2m f F M (M 2 ) , 
WEDM = dj = ie F^{M 2 Z ) . 

The anomalous weak-magnetic dipole moment (AWMDM) is the analogue of the anoma- 
lous magnetic dipole moment (AMDM) aj = (gj — 2)/2. The Fm (Fe) form factors are 
the coefficients of the chirality-flipping term of the CP-conserving (CP- violating) effective 
Lagrangian describing Z-fermion couplings. Therefore, they are expected to get contri- 
butions proportional to some positive power of the mass of the fermions involved. This 
allows the construction of observables which can be probed experimentally most suitably 
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by heavy fermions. Hence, for on-shell Z bosons, where the dipole form factors are gauge 
independent, the b quark and r lepton are the most promising candidates. 

The contribution to the AWMDM, for the r lepton and for heavy quarks, at the 
one-loop level has recently been calculated by Bernabeu et al. in the SM [|], || and in 
two Higgs doublet models . In this letter, we extend these calculations to the MSSM 
and investigate the maximum size of the non-standard contributions. We classify all the 
triangle diagrams in six topologies (classes) or generic diagrams (Fig. [I]) as in Ref 
Working in the 't Hooft-Feynman gauge all the would-be-Goldstone bosons are included. 
The global result, adding all the diagrams, is gauge independent. Every class of diagrams 
is calculated analytically for general couplings and expressed in terms of standard one-loop 
integrals C |§ and vertex coefficients A as follows: 

• Class la : C = C(—q,q,m,m,M), 

Fm{s) = ^-{4m f X v {-2C+ + 3C 1 -C ]+4mX v {-2C 1 +C ]}, (2) 

with X v = V{V' 2 + A' 2 ) ± 2AVA', X' v = V{V' 2 - A' 2 ). 

• Class lb : C = C(—q,q,mf,mf,0) 

*m(*) = C F ^{V[4m f (~2C+ + 3C 1 -C )+4m(-2C 1 +C )]}. (3) 

47T 

The colour factor Cp — |. 

• Class II : C = C(-q,q,M,M,m), 

F l A{s) = —{2m f \ v [ACt + C+}-Qm\' v Ct}, (4) 

with A y = G{V 2 + A 2 ), X' v = G(V 2 - A 2 ). 

• Class HI : C = C(—q,q,mx,m 2 ,M), 

*ff W = ^{-2 m/ Af [2CJ - C+] + (ft, + W + (ft, + ft 2 )A^ Cl - 

+ (wi - m 2 )\" A C+ + (mi - m 2 )X' v C^} , (5) 

with A^ = V(SS'-PP')±A(SP'-PS'), X'v = V{SS' + PP'), X'X = A{SP' + PS'). 

• Class IV : C = C(-q,q,M 1 ,M 2 ,rh), 

Fm(s) = ^{2m f X v \2C+ -C+}+mX' v [2Ci -C Q }} , (6) 

with Ay = G{SS' - PP'), X' v = G{SS' + PP'). 
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• Class V : C = C(-q,q,M 1 ,M 2 ,m), 

F&s) = ^{AWC^ + CT]}, (7) 

with A y = G(SV + PA). 

• Class VI : C = C(-q,q,M 1 ,M 2 ,m), 

F£(s) = 

with A y = G(SV-PA). 

We use the conventions of Refs. H for the couplings, masses and vertex coefficients.0 We 
employ the compact notation ay[ijk] for the contribution to the AWMDM of particles 
i,j, k running in the loop as labelled in Fig. C.l of Ref. ||. Notice that all the masses 
appearing explicitly as factors of the C integrals in the expressions are fermion masses, 
consistently with the chirality flipping character of the dipole moments. 

The WMDM of the r lepton 

In the SM there are 14 diagrams at one-loop, in the 't Hooft-Feynman gauge, contributing 
to the t-Z boson coupling. We are in agreement^ with Ref. with the total value for 
M H o = 2M Z , m T = 1.777 GeV, M z = 91.19 GeV, s 2 w = 0.232 and a = 1/128: 

af [SM] = (2.10 + 0.61 i) x KT 6 . 

The diagrams with Higgs bosons have only a small impact, changing the final result || 
by less than 1% for 1 < M H o /M z < 3. 

The MSSM contributions introduce new parameters, the values of which are not known 
but are constrained by present experiments. To reduce the number of free parameters, we 
assume the usual GUT constraint for the soft-breaking mass terms in the gaugino sector 
and a common squark mass parameter as well as a common slepton mass parameter. For 
given tan /3, we are left with the following free parameters: the gaugino mass parameter 
for the SU(2) sector M, the Higgs-higgsino mass parameter /i, the common slepton soft- 
breaking mass parameter mj, the common squark soft-breaking mass parameter rriq, the 
trilinear soft-breaking term A T and the mass of the pseudoscalar Higgs boson Ma- We 
consider two typical scenarios: low and high tan/?, respectively tan/? =1.6 and 50. The 
genuine MSSM diagrams to be included are: diagrams with MSSM Higgs bosons; diagrams 
with charginos and scalar neutrinos; and diagrams with neutralinos and f sleptons. 

1 This leads to the different signs of the tree level part of Eq. ([j]) when compared to Eq. (3.14) of 
Ref. §. 

2 The global sign is due to our different conventions. 



a 
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(8) 
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MSSM Higgs contribution to a™ 

The diagrams for Ztt involving Higgs bosons are: 
class III: [rrh], [ttH], [ttA], [ttG°], [wff-], [vvG~\, 

class IV: [Ahr], [G°hr], [AHr] , [G°Hr] , [hAr], [hG°r], [HAr] , [HG°r] , [H-H~v], 
[H-G-v], [G-H-p], [G-G-u], 

class V: [Zhv], [ZHu], [WG~v\, 

class VI: [hZv], [HZv], [G~Wv\, 

where G°, G ± are the would-be-Goldstone bosons and h, H, A and are the physical 
MSSM Higgs bosons. Actually, not all of them give a contribution to the AWMDM, in 
particular, the diagrams of class IV with neutral Higgs bosons vanish as in the case of the 
SM. The masses of the Higgs bosons are fixed by the mass of the pseudoscalar, Ma, tan/3 
and the common squark mass parameter rriq, which are the only free parameters involved 
in this calculation. A value rriq = 250 GeV will be assumed in the following. Low values of 
Ma yield the largest contribution and for Ma <^ 300 GeV the contributions become Ma- 
independent. The class III diagrams provide the only contribution to Im(a^), assuming 



the present lower bound for the masses of the Higgs bosons Due to the fact that the 
ratio of vector couplings v T /v u is small, the diagrams with charged Higgs bosons give the 
main contribution, which is of the order a/47r(m T /Mz) 4: tan 2 /3 ~ 10~ 10 (10~ 7 ) for the 
low (high) tan {3 scenario. For example, one gets 

Im(af [Higgs]) = -0.0001 (0.06) x 10~ 6 for M A = 100 GeV and low (high) tan/3 . 

For Re(ajf), the diagrams of classes III and IV are proportional to {m T /Mz) A tan 2 j3 and 
the diagrams of classes V and VI are proportional to (m T /Mz) 2 - Actually, 

^ , WrTT 1N f —0.3 (—0.4) x 10~ 6 for M A = 100 GeV „ , . , , 

Re a: Higgs = < „ tor low high tanrt . 

t 66 | _ Q g ( _ Q 3) x 1Q _ 6 for Ma > ^ V 6 ) 

Chargino and scalar neutrino contribution to a 1 ^ 

There are two diagrams involving charginos and scalar neutrinos: 
class III: [x7xjv], 
class IV: [uuXk]- 
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Charginos [12] and scalar neutrinos [13] do not occur in Z decays and hence their con- 
tribution to the AWMDM is real, in the region of their masses not ruled out by the 
experiments. The free parameters involved here are tan/3, M, /i and mj. The contribu- 
tion becomes smaller increasing mj (consistently with decoupling [JI0|). In the large M 
and \fi\ region the charginos also decouple. The contribution is enhanced by increasing 
tan /3. In the low tan (3 scenario the contribution is of order 10~ 7 . The value for mj = 250 
GeV, M = 200 GeV and |//| = 200 GeV is 

n n f -0.2 x 10~ 6 < 0) c , 

Re(a_ charginos ) = < . for low tanp . 

L J | +0.2 x 10~ 6 (fi > 0) 

In the high tan /3 scenario the contribution is of order 10 -5 . We find: 

n , Wn n f -7.0 x 10~ 6 (u<0) r , , 

Re[al charginos ) = < „ for high tanp 

6 J; \ +7.0 x 10- 6 {/I > 0) S ^ 

and the same values for M and \fi\ as above. This enhancement for higher tan /3 is expected 
in analogy to the AMDM of the muon. The contours in the M — fi plane are similar to 
the ones for the b quark (Fig. 0) suppressed by factors between m T /m b and (m T /mb) 2 [see 
Eqs. (H) and 



Neutralino and f slepton contribution to 

There are two diagrams involving neutralinos and f sleptons: 
class III: [XiXjTk], 
class IV: [nfjxl]- 

The f slepton masses are bounded by present experiments ]14[ to be heavier than half 



the mass of the Z, but the neutralinos are less stringently bounded and could be much 



lighter [12], allowing for the possibility of a contribution to the imaginary part of the 
AWMDM through diagrams of class III. Now, we have tan/3, M, /i, m z ~, as well as the 
trilinear soft-breaking mass term for the f sleptons, A T , as free parameters. We choose 
different values for A T in the range ±/itan/3 to estimate the sensitivity of the AWMDM 
to this soft-breaking term and find no sizeable deviation. However, for some values of 
A T and [i the mass of the f sleptons can become unphysical. For simplicity, we take the 
value of A T that makes m T LR = A T — /itan/3 = 0. Just below the neutralino threshold the 
imaginary part is 

| Im(a^ [neutralinos]) | ~ 10~ 9 (10~ 7 ) for low (high) tan/3 
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and it is negligible far below threshold. Concerning the real part, also the value is max- 
imum near the neutralino pair threshold. Like in the chargino-scalar neutrino contribu- 
tions, the result is expected to be enhanced by tan (3. We find 

Re(af [neutralinos]) = -0.02(-0.4) x KT 6 for low (high) tan/3 

and M = fi = 200 GeV. 

The WMDM of the b quark 

In the electroweak SM there are 14 diagrams, in the 't Hooft-Feynman gauge, contributing 
to the b-Z boson coupling at the one-loop level, plus one more (gluon exchange) if QCD is 
also included. Our calculation reproduces the result of Ref. j|. Taking as input m h = 4.5 
GeV, m t = 175 GeV, M z = 91.19 GeV, s 2 w = 0.232, a = 1/128 and a s = 0.118, 
the pure electroweak contribution is of(EW) = [(1.1; 2.0; 2.4) - 0.2 i] x 10" 6 , for 
M H o = M z , ^M z , SM Z , but the whole value is dominated by the QCD contribution to 
give 

ag [EW + QCD] = (-2.98 + 1.56 i) x 10" 4 . 

Only the third generation of quarks is considered and the element V t b of the CKM matrix 
is set to 1. 

In the MSSM one has to include the following genuine supersymmetric contributions: 
diagrams with MSSM Higgs bosons; diagrams with charginos and t squarks; diagrams 
with neutralinos and b squarks; and diagrams with gluinos and b squarks. In this case, 
we are left with the following free parameters: M, /i, rriq, M^, the trilinear soft-breaking 
terms and A t and a mass for the gluinos m§. 

MSSM Higgs contribution to ag 

The diagrams for Zbb involving Higgs bosons are 
class III: [bbh], [bbH], [bbA], [bbG°], [UH~], [ttG~], 

class IV: [Ahb], [G°hb], [AHb] , [G°Hb], [hAb], [hG°b], [HAb], [HG°b] , [H~H~t], [H~G~t], 
[G~H~t], [G~G~t], 

class V: [Zht], [ZHt], [WG~t], 

class VI: [hZt], [HZt], [G~Wt]. 
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As before, the diagrams of class IV with neutral Higgs bosons do not contribute. The 
diagrams of class III with neutral Higgs bosons yield Im(a^) ~ a/47r(mb/M^) 4 tan 2 8 ~ 
10 -8 ( 10 -6) f or i ow (^gh) tan/?. Taking m~ q = 250 GeV we get 

T , w rTy n f °- 02 ( 3 - 7 ) x 10-6 for m a = 100 GeV , , 

lmial Higgs = < „ and low high tanp . 

V b 1 U \ 0.01 (0.08) x 10~ 6 for M A = 300 GeV V ; 

Concerning Re(a^), there are contributions proportional to > ^ m 2 ~ ) an< ^ 



( ) tan 2 3. The sum is thus not very sensitive to tan 8, 

\Mz) 



,, rT , | -3.8 (+0.8) x 10- 6 for M A = 100 GeV . . 

Re ar Higgs = < . and low high tanp . 

1 1 1 -1.8 (-0.9) x 10- 6 for M A = 300 GeV V ; 

Chargino and t squark contribution to 

There are two diagrams involving charginos and t squarks: 
class III: \XiXjik], 
class IV: [i{tjX^]. 



As in the r case, assuming the present bounds on chargino and squark masses |12| , 15[ 
lm(a^ [charginos] )=0. The free parameters involved are tan/?, M, /i, rriq and A t . Varying 
the value for A t in the range i/icot/3 shows that there is only a small effect in the 
AWMDM. We thus choose A t in a way that the off-diagonal entry of the t squarks mixing 
mass matrix vr^ LR — A t — \i cot 3 = 0. Increasing rriq the contribution becomes smaller. 
Fig. [2] shows a contour plot in the M — fi plane for rriq = 250 GeV. In the large M and 
region, the charginos decouple, as expected from the decoupling theorem [|Tl]]. The 
contour lines for the masses of the lightest chargino and lightest neutralino in the M — \i 
plane are also depicted. The contributions are enhanced by increasing tan f3. In the low 
tan/3 scenario the contribution is of order 10~ 6 (Fig. |]a). For M = =200 GeV we get 

"r / wr l • n / -1.1 x 10~ 6 (^ < 0) 
Heidi charginos ) = < a for low tanp . 

1 b 1 S J; \ 1.0 x 10- 6 (/i > 0) ^ 

In the high tan/3 scenario (Fig. 0b) the contribution is enhanced to become of order 10~ 5 . 
For the same values of M and fi as above, 

1 1 ■ , i -32.2 x 10- 6 (p<0) 

Heidi charginos I ~ < „ for high tanp . 

" [ ^ f 33.6 x 10~ 6 (im > 0) S ^ 
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Neutralino and b squark contribution to 

There are two diagrams involving neutralinos and h squarks: 
class III: [xUfh], 
class IV: \bibjxV\- 

Now, we have tan/3, M, well as the trilinear soft-breaking mass term for the 

b squarks, Aj,, as free parameters. There is no sizeable dependence on A^. As in the 



t case, due to the heavy mass of the b squarks fll5| , an imaginary part can only arise 
through diagrams of class III. The region of non vanishing Im(a^) is below the neutralino 
threshold and the value is 

| Im(aj^ [neutralinos] )| ~ 0.1(1) x 1(T 6 for low (high) tan/3 , 

near threshold. Concerning the real part, also the value is maximum near the neutralino 
pair threshold. The result is proportional to tan (3 as expected: 

Re {af [neutralinos]) = -0.2 (-10.5) x 1(T 6 for low (high) tan/3 

and M = /i = 200 GeV. 

Gluino contribution to 

There is one diagram with a gluino and two b squarks running in the loop: 
class IV: [bibjg]. 



The present bounds on b squark masses |15[ only allow a real contribution to the AWMDM. 
The free parameters are tan/3, rrig, m g , Af, and \i. These last two parameters affect the 
result only through the off-diagonal term of the b squark mass matrix m b LR = Af, — fi tan /3. 
There still remains a slight dependence on tan /3 in the diagonal terms of the b squark 
mass matrix and so we keep the distinction between low and high tan/? scenarios. The 
mixing in the b sector is determined by m b LR and intervenes in the contribution due to 
chirality flipping in the gluino internal line (the contribution proportional to rrig). Thus, 
the contribution to the AWMDM is enhanced by the largest values of m b LR compatible 
with a physical and experimentally not excluded mass for the lightest b squark. Assuming 
a typical scale for m b LR of order /ztan/? leads to values of the order 100 GeV (10 TeV) for 
the low (high) tan /3 scenarios. 

The low tan/3 region is displayed in Fig. [3]a, for rrig = 250 GeV. The mixing for the 
plotted values of m b LR do not affect significantly the mass of the lightest b squark (always 
above 200 GeV). For zero gluino mass, only the term proportional to the mass of the b 
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quark provides a contribution. As we increase the gluino mass, the term proportional to 
nig dominates, especially for large m b LR , being again suppressed at high m g due to the 
gluino decoupling. The typical value is 

Re(a^[gluinos]) = -1.1 x 1(T 6 for low tan/3 

and m g = m b LR = 200 GeV. For high tan/5 (Fig. ^|b), the behaviour is analogous but 
larger values can be obtained. For example, choosing m g = 200 GeV and m b LR = 7 TeV, 
we obtain 

Re(a^[gluinos]) = -18.9 x 10" 6 for high tan/3 . 

Conclusions 

From the previous analysis we conclude that the imaginary part of the AWMDM of the r 
lepton is provided only by the diagrams with Higgs bosons, being at most of order 10~ 7 
for the high tan/3 scenario, an order of magnitude below the SM contribution. For the 
real part, the charginos dominate with a contribution of order 10~ 5 (10 -6 ) in the high 
(low) tan/3 scenario at the level of the SM contribution or even larger. The neutralino 
contribution is of opposite sign, but one order of magnitude smaller. Also the Higgs 
contribution to the real part is negligible. 

The most important MSSM contribution to the is provided by charginos and 
gluinos in the high tan (3 scenario. For high tan /3, the neutralino and chargino contribu- 
tions have opposite sign, but the former is typically a factor five below. On the other 
hand, the chargino contribution has the sign of /i whereas the gluinos only contribute with 
negative sign. Therefore the total contribution to the real part is maximal for negative \x 
and can reach the value of 

Re(af [MSSM]) ~ -60 x 10~ 6 for high tan/3 . 

These values are one order of magnitude higher than the pure electroweak SM contri- 
bution, but still a factor five below the standard QCD contribution. The Higgs boson 
diagrams are negligible, and they are only important for the imaginary part, contributing 
the same amount as the neutralinos below threshold, at most 

|Im(a^[MSSM])| ~ 10~ 6 for high tan/3 . 

This value is of the same order as the electroweak SM result. In the low tan/3 scenario, 
all the MSSM contributions to the real part are of order 10~ 6 ; for the imaginary part 
they are one order of magnitude smaller. They both are comparable in size to the pure 
electroweak SM contribution. 
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Class la: Class lb: 




Figure 1: The one-loop Zff diagrams with general couplings. 
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Figure 2: Re(a^[charginos]) for low (a) and high (b) tan j3 in units of 10 -6 in the plane M — (i. 
The contour solid-lines correspond to the light chargino masses m^±=80 and 200 GeV and the 
contour dashed-lines correspond to the lightest neutralino masses m^o=15 and 100 GeV. The 
common squark mass parameter is fixed to rriq = 250 GeV and m t LR = 0. 
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Figure 3: Re(a^ [gluinos]) low (a) and high (b) tan/3 in units of 10 in the plane m LR — rrig. 
The common squark mass parameter is fixed to rriq = 250 GeV. The region above the upper 
horizontal line is unphysical (m? < 0). 
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